Mesoscopic two-phase model for describing apparent slip in micro- channel flows 
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The phenomenon of apparent slip in micro-channel flows is analyzed by means of a two-phase 
mesoscopic lattice Boltzmann model including non-ideal fiuid-fluid and fluid-wall interactins. The 
weakly-inhomogeneous limit of this model is solved analytically. The present mesoscopic approach 
permits to access much larger scales than molecular dynamics, and comparable with those attained 
by continuum methods. However, at variance with the continuum approach, the existence of a gas 
layer near the wall does not need to be postulated a priori, but emerges naturally from the underlying 
non-ideal mesoscopic dynamics. It is therefore argued that a mesoscopic Lattice Boltzmann approach 
with non-ideal fluid-fluid and fluid-wall interactions might achieve an optimal compromise between 
physical realism and computational efficiency for the study of channel micro-flows. 
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The microscopic physics underlying fluid/solid inter- 
actions is fairly rich and complex, for it depends on spe- 
cific details of molecular interactions as well as on the 
micro-geometrical details of the boundary. However, on 
a macroscopic scale, these details can often be safely ig- 
nored by assuming that the net effect of surface interac- 
tions is simply to prevent any relative motion between the 
solid walls and the fluid elements next to them. This is 
the so-called "no-slip" boundary condition, which forms 
the basis of mathematical treatments of bounded flows 
as continuum media . No-slip boundary conditions are 
extremely successful in describing a huge class of viscous 
flows. Yet, the evidence is that certain classes of viscous 
flows do slip on the wall. Recent advances in microfluidics 
experiments 0) 01 j as well as numerical investigations 
0,1^0,13,13, have identified the conditions which seem to 
underlie the validity of the no-slip assumption. Namely: 
(i) single-phase flow; (ii) wetted surfaces and (iii) low lev- 
els of shear rates. Under such conditions, careful exper- 
iments have shown that fluid comes to rest within a few 
molecular diameters from the surface 

iliailllll Con- 
ditions (i-iii) are not exhaustive, though. For instance, 
partial slips of simple (Newtonian) flows, such as alkanes 
and water, is pred icted by an increasing number of ex- 
periments [H 111 El m and simulations H H E S Q 
(see [13 for a review on experiments and numerics) . Un- 
der this state of affairs, there appears to be a great need 
to provide a convincing, and possibly general, theoreti- 
cal picture for the onset of slip motion. Among others, 
an increasingly popular explanation is that the flowing 
fluid would develop a lighter (less dense) phase and dy- 
namically segre gate it in the form of a thin film sticking 
to the wall [13, uM ■ This thin film would then provide 
a "gliding" surface for the bulk fiuid which would slip 
on it without ever coming in contact with the solid wall. 
This gives rise to the so-called apparent slip phenomenon, 
that is, the extrapolated bulk fiow speed would vanish 
far-out away from the wall, even though the actual fiow 



speed in the film does vanish exactly at the wall location. 
This film-picture is very appealing, but still in great need 
of theoretical clarification. In particular, the underlying 
mechanisms of film formation are still under question: 
are they generic or detail-driven? 

In this paper we shall propose that film formation is a 
generic phenomenon, which can be captured by a one- 
parameter mesoscopic approach, lying in-between the mi- 
croscopic (atomistic) and macroscopic (continuum) lev- 
els. The mesoscopic approach is based on a minimal 
(lattice) Bohzmann equation (LBE) [20II2II2II includ- 
ing non-ideal interactions [23 . 12I iMlM l27l l28l| . which 
can drive dynamic phase transitions. The only free pa- 
rameter in the LBE is the strength of these non-ideal 
(potential energy) interactions. Hopefully, the present 
mesoscopic approach provides an optimal compromise 
between the need of including complex physics (phase- 
transition) not easily captured by a continuum approach, 
and the need of accessing experimentally relevant space- 
time scales which are out of reach to microscopic Molec- 
ular Dynamics (MD) simulations j^, iSi iS] • I^i particu- 
lar, at variance with the macroscopic approach, the gas 
film does not need to be postulated a-priori, but emerges 
dynamically from the underlying mesoscopic description, 
by progressive switching of potential interactions. One 
major advantage of this formulation is that it allows to 
develop a simple and straightforward analytical interpre- 
tation of the results as well as of the effective slip length 
arising in the flow. This interpretation is based on the 
macroscopic limit of the model which can be achieved by 
a standard Chapman-Enskog expansion. 
The lattice Boltzmann model used in this paper to de- 
scribe multiple phases has been developed in |23|. Since 
this model is well documented in the literature, here we 
shall provide only the basic facts behind it. We recall that 
the model is a minimal discrete version of the Boltzmann 
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equation, and reads as follows: 

fi{x + ci,t+l)- Mx,t) = (^fiix,t) ~ fl^'^x.t)) 

(1) 

where fi{x,t) is the probability density function associ- 
ated to a mesoscopic velocity c; and where r is a mean 
collision time and fj^'^'^\x, t) the equilibrium distribution 
that corresponds to the Maxwellian distribution in the 
fully continuum limit. The bulk interparticle interaction 
is proportional to a free parameter, (yb, entering the bal- 
ance equation for the momentum change: 

^^ = F = gfc^«;zVl,[p(a;)]vl/[p(a;-KQ)]cz (2) 

I 

being wi the equilibrium weights and ^' the potential 
function which describes the fluid-fluid interaction trig- 
gered by density variation. By Taylor expanding eq.Q 
one recovers, in the hydrodynamical limit, the equation 
of motion for a non-ideal fluid with equation of state 
P — c^l^p— ^Qb'i'^ip)), Cs being the sound speed velocity. 
With the choice 
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FIG. 1: Rarefaction effects in the full-interaction case 
QwiGb 7^ 0. Density profiles are plotted as a function of 
the distance from the wall, normalized to the channels height 
(y/Ly). The wall interactions have been fixed assuming 
Qw ~ 0.03 and ^ = 2. The following values of Qt are con- 
sidered: Gt = 1.5 (+), Gb = 2.5 (x), Gb = 3.5 (★), Gb = 3.9 
(□), Gb ~ 3.98 (o). We remind that the bulk phase transition 
is set at Gw = 4. In all simulations we choose Ly — 80 grid 
points for the height of the channel. The volume averaged 
Knudsen is Kn ~ 10""^, this would correspond to a channel 
of a few fim for liquid water. 



^{p) = 1 - exp(-p/po) 

with po = 1 a reference density, the model supports phase 
transitions whenever the control parameter exceeds the 
critical threshold Qt, > 0§. In our case, tyj = 4 for an 
averaged density (p) = log(2). 

We consider Qh as an external control parameter, with no 
need of responding to a self-consistent temperature dy- 
namics. It has been pointed out [s^l that the SC model 
is affected by spurious currents near the interface due to 
lack of conservation of local momentum. This criticism, 
however, rests on an ambiguous interpretation of the fluid 
velocity in the presence of intermolecular interactions. In 
fact, spurious currents can be shown to disappear com- 
pletely once the instantaneous pre and post-coUisional 
currents are replaced by a time-average over a collisional 
time. This averaged quantity is readily checked to fulfill 
both continuity and momentum conservation equations 
without leading to any spurious current |3l[ . Let us now 
consider the main result of this letter, namely the critical 
interplay between the bulk physics and the presence of 
wall effects. In fact, in order to make contact with exper- 
iments and MD simulations, it is important to include 
fluid-wall interactions, and notably a parametric form 
of mesoscopic interactions capable of mimicking wetta- 
bility properties as described by contact angles between 
droplets and the solid wall [33 • This effect is achieved 
by assuming that the interaction with the wall is repre- 
sented as an external force F^ ,, normal to the wall and 
decaying exponentially (2^113, i.e. 

F,,{x) = g^p{x)e-\^—\/i (3) 

where x^ is a vector running along the wall location and 
^ the typical length-scale of the fluid-wall interaction. 



Equation J^Jl has been previously used in literature by 
using a slightly different LBE scheme to show how the 
wetting angle depends on the ratio Qw/Qb in presence of 
phase coexistence between vapor and liquid |28|. Here 
we want to study the opposite situation, i.e. the effects 
of Qw when the thermodynamically stable bulk physics 
is governed by a single phase. The main result is that 
the presence of the wall may trigger a local phase co- 
existence inducing the formation of a less dense phase 
in the vicinity of the walls and an apparent slip of the 
bulk fluid velocity profile extrapolated at the wall loca- 
tion. Equations H1I3|) have been numerically solved for 
different values of the parameters t/b, Qw and ^ in a two 
dimensional channel with periodic boundary conditions 
in the stream-wise x direction, being y — Q and y = Ly 
the wall positions. The sign of Qw is such to give a re- 
pulsive force for the liquid particles at the wall. 

The flow is driven by a constant pressure gradient in 
the X direction Fi — Si^xdxPo- No-slip boundary condi- 
tions are used at the wall and for small Knudsen numbers, 
i.e. in the large scale limit, the numerical solutions have 
been checked against its weakly-inhomogeneous macro- 
scopic hydrodynamic limit, namely: 

dtp + d,iu,p) = (4) 
p [dtu^ + {ujdj)ui] = ~diP + vdj{pdiU.j + pdjUi) + Fi 
P^clp-V,ff{p) 

where subscripts i, j run over the two spatial dimensions. 
Above we have v = (?s{t — 1/2) and P is the total pres- 
sure consisting of an ideal-gas contribution, c^p, plus the 
so-called excess pressure, K//, due to potential-energy 
interactions. The expression of Ve// in terms of both Qi, 
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and reads: 



Veffip) ^ ^Gbil 



dsp{s) exp(— s/^). 



Let us notice that the continuum equation Q naturally 
predicts the increase of the mass flow rate in presence 
of a density profile which becomes more and more rar- 
efied by approaching the wall 1^. Indeed, under station- 
ary conditions, the continuity equation in Q reduces to 
dy{pUy) = 0, which, because of the boundary conditions, 
implies puy = 0, i.e. Uy = Q everywhere. Thus, in a ho- 
mogeneous channel along the stream-wise direction, the 
velocity satisfies the equation 



vdy{pdyUx) = -dxPa- 



(5) 



In the new variable, y' = y — H, where H — Ly/2, we 
may express the solution of jSJ as: 



Ux{y') 



' vp{s) 



ds. 



(6) 



Using ^ and assuming that density variations are con- 
centrated in a smaller layer of thickness 5 near the wall, 
we can estimate the mass flow rate Qeff for small 5 as: 



3 A/9„, 5 



(7) 



where Qpois corresponds to the Poiseuille rate 
"^dxPoH^ /iv valid for incompressible flows with no-slip 
boundary conditions. In equation |(7J), the quantity Ap^ 
is defined as the difference between p computed in the 
center of the channel and pw computed at the wall. The 
effective slip length is then usually defined in terms of 
the increment in the mass flow rate jlTjj: 



A, 



(8) 



This is the best one can obtain by using a purely con- 
tinuum approach. The added value of the mesoscopic 
approach here proposed consists in the possibility to di- 
rectly compute the density profile, and its dependency on 
the underlying wall-fluid and fluid-fluid physics. To this 
purpose, we consider the momentum balance equation in 
for the direction normal to the wall, i = y. Since 
Uy — 0, we simply obtain dyP = 0, i.e. 

cldyp - 206(1 - e-P)e-PdyP - g.^pe~y/^ = 0. (9) 

Let us first study the effects of the wall in © by set- 
ting Qb = 0. One can easily obtain \og{p{y)/p^) = 
^^(1 — exp(— j//^)), which enables us to estimate Ap^, — 

Pw{s^p{iGw/cl) — 1). Using (|SJ), we obtain for the effec- 
tive slip-length: 



[Gb = 0]. 



(10) 
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FIG. 2; Increase of the mass flow rate with the coupling 
strength Qt of fluid-fluid bulk interactions. Fixing the wall 
correlation function ^ = 2, we plot the mass flow rate {Qeff) 
normalized to its Poiseuille value {Qpois) as a function of C/t 
for different values of Qwaii- Gwaii = 0.0 (o), g^^aii = 0.04 (□), 
Gwaii = 0.08 (x). Inset: same as the main figure for Ap^/p. 



We now turn our attention to the non trivial interference 
between bulk and wall physics whenever Gb > 0. Defining 
the bulk pressure as: Pb — c^sP ~ \Gb{^ — exp(— p))^, we 
can rewrite equation (j^l to highlight its physical content 
as follows: 



log 



(piyl 

V Pw 



^GA^-e'y/i)/dPb/dp (11) 



where the bulk effects appear only through the following 
term: 



1 



dPb ^ 

dp ~ \og{p{y)lpu 



y dPb dp 
dp p ' 



(12) 



Equation (fTT|l highlights two results. First, the effect of 
the bulk can always be interpreted as a renormalization 
of the wall-fluid interaction by 



G^l 



dPb 
dp 



(13) 



Second, as it is evident from (|13|) . one must notice that 
near the bulk critical point where dPb/dp — > 0, the renor- 
malizing effect can become unusually great. In other 
words, the presence of the wall may locally push the sys- 
tem toward a phase transition even if the bulk physics it 
is far from the transition point. As a result, the effective 
slip length in presence of both wall and bulk non-ideal 
interactions can be estimated as: 



A, ~Cexp(C0«) 



(14) 



In Fig. n we show p{y) for different values of Gb and 
Gw = 0.03, ^ = 2 as obtained by numerically in- 
tegrating equations H1I3|) . The numerical simulations 
have been carried out by keeping fixed the value of 
(p) = Y~ Iq^ p{s)ds = log(2). As one can see, while 
Gb Gc = 4, the density difference Ap^; between the 
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FIG. 3: Momentum profile as a function of tiie cliannel lieiglit. 
We plot the momentum profile i^pUx) normalized to its center 
channel value ((pua;)c) as a function of the distance from the 
wall normalized to the channel height (Ly). The results 
of numerical simulations (□) with Qh = 3.5, = 0.08 and 
^ = 2 are compared with the analytical estimate (continuous 
line) obtained solving equations @ and 0. To highlight the 
rarefaction effect, the parabolic fit in the center channel region 
(dotted line) is also plotted. Inset: estimate of the apparent 
slip length in the channel obtained the same parabolic fit as 
in the main figure. 



center of the channel and the wall increases, as predicted 
by equation @. Consequently, the mass flow rate in- 
creases as shown in Fig. |21 Let us notice in the same 
figure that also with = 0, the wall initiates a small 
rarefaction effects due to the fact that fluid particles close 
to the boundary are attracted only by particles in the 
bulk of the channel. What we showed here is that the 
combined actions of and Q\, may strongly in- 

crease the formation of this less dense region in the prox- 



imity of the surface. For a quantitative check, we have 
numerically integrated equations Q and © for a given 
value (p) = log(2). The analytical estimate for pux is 
compared with the numerical results in Fig. |21 This is 
a stringent test for our analytical interpretation. The 
result is that the analytical estimate is able to capture 
the deviations from a pure parabolic profile at approach- 
ing the wall region, where rarefaction effects are present. 
The crucial point in our analysis is that, even for very 
small large apparent slip can occur in the channel if 
Qh is close to its critical value, i.e. the limit Qw — > and 
Qh Qb do not commute. For example, let us consider 
the case when t/^ ~ e <C 1, ^ ~ e and Qt ^ Q^ — e"^; 
we obtain e'^ and therefore, equation predicts 

that As ~ 0(1) for e ^ 0. The wall effect, parametrized 
by Qj^ and ^, can act as a catalyzer in producing large 
apparent slip. Most of the results shown in Figs. (Q) and 
(0) are in close agreement with the MD numerical sim- 
ulations 0,0, HQ- Our analysis points out that, close 
to the wall, one can observe a "local phase transition" 
triggered by the presence of the wall itself. In summary, 
we have shown that a suitable form of the Lattice Boltz- 
mann Equation can be proposed in order to simulate ap- 
parent slip in microchannel. Slip boundary conditions 
arise spontaneously because, close to the wall, a "gas" 
layer is formed. If the system is close to a state where 
coexistence of different phases (liquid and gas) are ther- 
modynamically achievable, then, macroscopic slip effects 
can result. We have shown that for large scale separation, 
the model reduces to a continuum set of hydrodynamical 
equations which explains the qualitative and quantita- 
tive behavior of the mass flow rate in terms of the model 
parameters, i.e. Qb and Qw 
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